Nonlinear Schrodinger equation for a superfluid Fermi gas from BCS to Bose crossover 
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We introduce a quasi-analytic nonlinear Schrodinger equation with beyond mean-field corrections 
to describe the dynamics of a zero-temperature dilute superfluid Fermi gas in the crossover from 
the weak-coupling Bardeen-Cooper-Schrieffer (BCS) regime, where kp\a\ -C 1 with a the s-wave 
scattering length and k F the Fermi momentum, through the unitarity limit, k F a — > ±00, to the 
Bose regime where fea > 0. The energy of our model is parametrized using the known asymptotic 
behavior in the BCS, Bose, and the unitarity limits and is in excellent agreement with accurate Green 
function Monte Carlo calculations. The model generates good results for frequencies of collective 
breathing oscillation of a trapped Fermi superfluid. 



PACS numbers: 71.10.Ay, 03.75.Ss,67.85.Lm, 05.30.Fk 



The crossover from a Bardeen-Cooper-Schrieffer (BCS) 
Fermi superfluid at zero temperature for weak coupling 
to a Bose condensate of dimers [l[ has been an intense 
area of research (both experimental @, Q and theoretical 
0, B H 0, B S El) after the realization of a BCS to 
Bose crossover (BBC) in a trapped dilute Fermi super- 
fluid near a Feshbach resonance. This allows to change 
the system from the BCS regime, with small negative 
a (/cfM *C 1), through the unitarity regime of diver- 
gent a {kpa - * ioo), to the Bose regime of dimers with 
positive a, with kp the Fermi momentum. The limiting 
behavior of the system in the weak-coupling (fcj?|a| <C 1) 
PH, E3, Q 51 and unitarity limits for both posi- 
tive and negative a are well known 

[HG1G1- Accu- 
rate information about the BBC dynamics has recently 
been available from numerical fixed-node Green's func- 
tion Monte Carlo (GFMC) calculations [1, 0] at zero 
temperature. There is also mean-field BCS (MFBCS) 
calculation of the same [6[. 

A quasi-analytical model for the BBC problem, gen- 
erating a the Ginzburg-Landau (GL) type equation for 
fermions for a < and a Gross-Pitaevskii (GP) type 
equation for dimers for a > 0, both including beyond 
mean-field effects, so as to be valid in the unitarity limit 
of divergent a, should be useful. We propose such a model 
in three dimensions, called the BBC model (BBM), for 
the crossover problem of a dilute trapped Fermi super- 
fluid at zero temperature using the known theoretical so- 
lution in the weak-coupling and unitarity limits in both 
the BCS and Bose regimes. 

We consider a dilute superfluid Fermi gas of N spin- 
half atoms of mass to, atomic scattering length a, and 
density n with singlet pairing due to an attractive atomic 
interaction. We assume, for a dilute gas, that the results 
are universal determined solely by the scattering length a 
and independent of the detail of the interaction potential 
and its range. We present an analytical model for energy 
and bulk chemical potential in the entire crossover region 
and hence derive a nonlinear equation for the superfluid 
Fermi gas. We present results for radial and axial fre- 
quencies of collective oscillation in a cigar-shaped trap 
and compare with GFMC and MFBCS calculations as 
well as experimental data. 



At low densities, in the BCS regime (fc F |a| <^C 1), gaps 
are negligible 0, [lj] and the total energy E per particle 
of the superfluid Fermi gas is given by [Tj], UM 
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where E F = fr 2 k 2 F /(2m) = An 2 / 3 = A_ (3^2/3 ig the 
Fermi energy with kp = (?>'n 2 n) 1 / z 

In the Bose regime (a > 0), the paired fermions form a 
weakly repulsive Bose gas of dimers of mass to' = 2to and 
density n' — n/2 with a dimer-dimer scattering length 



a' = 0.6a, as predicted by Petrov et al. [1 
regime, the energy of the system is given by 
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with a' = 0.6a [17], where eb is the (positive) binding 
energy of the dimcr. The dynamics of the dimer bosons 
is governed by the energy (E/N + £s/2) discounted for 
dimcr binding. The lowest-order term in Eqs. |T|) and ([2|) 
leads to the standard mean-field GL equation for Fermi 
superfluid and GP equation for dimer bosons, respec- 
tively, with higher-order terms leading to beyond mean- 
field corrections due to interaction. 

The limiting energies ([T]) and @, valid for small |a|, 
drastically fail as one approaches the unitarity limit 
a — * ±00. For a dilute Fermi superfluid, the unitarity 
limit is supposed to be universal and the relevant en- 
ergy scale is the energy of the noninteracting Fermi gas 
3E p/5 [l4l . [H ] , so that the energy of the system is given 
by [HOI] E/N = 3Ept;/5 with £ a universal factor. By 
using Pade approximant Baker [13| estimated two val- 
ues for £ = 0.326 and 0.568, and Heiselberg [l2| obtained 
£ = 0.326. Later, from experimental results the following 
values for £ have been obtained f = 0.51 ± 0.04 UM , 0.7 
[U, 0.27j££g H3, 0.41 ±0.15 0, 0.46 ±0.05 \2$, and 
12 01 • The accurate theoretical estimates for the 
factor £ are £ = 0.44±0.01 and 0.437±0.009 (GFMC) 
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[3, and 0.59 (MFBCS) @]. Bulgac and Bertsch Qj| sug- 
gested the following approximate behavior of energy per 
particle near the (fopM) -1 — » limit 
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in both BCS (a < 0) and Bose (a > 0) regime with 
£ w v w 1. In the Bose regime (a > 0), one should 
replace § by § + ^f. 

We write a simple expression for energy E/N combin- 
ing the limiting behaviors dXJ) and §5§ in fermion (a < 0) 
and dimer boson (a > 0) regime. We suggest the follow- 
ing expression for fermion energy for a < 
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This expression has the constant 10/(97r) taken from 
limit ([1]) and the constant £ from limit ([3]). Expression 
((!]) reproduces the first two terms in Eq. (p} exactly for 
small \a\. It also reproduces the first term in Eq. ([3]) for 
large \a\, for the second term it yields ( = 0.9186 in place 
of the approximate value £ w 1 

For a > 0, we use the asymptotic behaviors of energy 
for large and small a, Eqs. ||3J| and ([3]), to propose the 
following expression for energy 



E e B 3 
JV 2~ 5 



:(*FO ; ) 



27V6tt 



= (fc F a') 5/2 



64(a7a)(fc F a , ) 3 / 2 64(fc F a') 5 / 2 



27£ 2 V6tt 5 



274V67T 5 



.(5) 



For large a', by construction, expression ([5]) satisfies Eq. 
© with £ = 1. For small a it reproduces the first term 
of expansion © exactly and the next terms closely. 

The only free parameter in model Eqs. ^ and (JSJ) 
(termed BBM) is the universal factor £, for which we use 
£ = 0.44 0,0, El- The k F a dependence of BBM © and 
([5]) are consistent with the behavior for small and large 
fcp|a|. Kim and Zubarev [9J considered two different [2/2] 
Pade approximants to parametrize energy in the Fermi 
superfluid and Bose dimer regimes. Manini and Salasnich 
!§j considered an energy function with arctan dependence 
on scattering length both in the Fermi superfluid and 
Bose dimer regimes. 

Next we plot in Fig. [TJ for a < 0, the energy from Eq. 
(U|) , in addition to those from limits (Q]) and © . We also 
plot the results of GFMC calculations of Refs. 0, 03 
and parametrizations of Refs. [1, Q. Next we plot, for 
a > 0, in Fig. [21 the results for energy from Eq. ([S]), the 
asymptotic limits © and (J3J, the GFMC results of Refs. 
0, [lj] and the parametrizations of Ref. [1, H| . 

From Figs. Q] and we realize that limits |T]) and 
([3j) , essentially, determine the energy for Uf\o\ < 1 and 
kp\a\ > 10. The correct energy over the entire crossover 
should be a smooth interpolation between these limits. 
This is what has been done to obtain the present results 
for a < and a > in very good agreement with the ac- 
curate GFMC calculations 0, Q3] and asymptotic limits. 




FIG. 1: (Color online) EHNEf) vs. fc F |a| for a < 0: BBM 
- Eq. ©; KZ from Ref. Q; MS from Ref. [|; limiting value 
(kFa)^ 1 — > —0 of Eq. ([3]); limiting value — > —0 of Eq. 
(TJ1); GFMC (1) [13; GFMC (2) 0. 
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FIG. 2: (Color online) E/(NE F ) +e B /(2E F ) vs. k F a for 
a > 0: BBM - Eq. ©; KZ from Ref. [|; MS from Ref. 
0; limiting value (fc^a)" 1 — ► +0 of Eq. @; limiting value 
k F a +0 of Eq. ©; GFMC (1) 0; GFMC (2) 0. 



To study the collective and bulk properties of the su- 
perfluid Fermi gas, we write convenient mean-field equa- 
tions. They are, essentially, the Schrodinger equation 
with a nonlinear term equal to the bulk chemical po- 
tential /Lt(n, a) = d£/dn, with £ = (E/N)n the energy 
density. A straightforward calculation with the leading 
two terms of energy of Eq. ([T]) yields the following bulk 
chemical potential for a — > — [4| 



(x(n,a) = An 1 l z + 2irh 2 an/m, 
to be used in the nonlinear Schrodinger equation 
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where U(r) is a trap, (r,t) are space and time variables, 
density n is related to order parameter ^ by n — \^\ 2 , 
such that J ndr — N, and m* ~ 2m is the effective mass 
of pairs. This equation can also be interpreted to be 
the Euler-Lagrange equation of a time-dependent den- 
sity functional theory [gL l25[| . Also, for a large number of 
atoms N, for a description of collective properties Eq. ([6|) 
is completely equivalent to the quantum hydrodynamic 
equations d . The first term on the right-hand side of Eq. 
([5]) is the standard nonlinearity in the three-dimensional 
GL equation of a Fermi superfluid (equal to Fermi en- 
ergy) with the last term appearing due to atomic inter- 
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action. This last term is reduced by a factor of 2 com- 
pared to the GP nonlinearity for bosons [Aixh 2 an / m) , as, 
in the present case, half of the atomic interactions, those 
between spin-parallel fermions, are inoperative. 

The bulk chemical potential p(n,a) corresponding to 
BBM energy gj) 



p{n,a) = An 2/3 - 
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is to be used in Eq. J7J). A simpler expression for bulk 
chemical potential can be obtained if we recall that, in the 
unitarity limit (a — » — oo), the energy is given by E/N = 
3Ep£/5 corresponding to a bulk chemical potential [l3[ 



p(n,a) = ?j 2 (37r 2 n) 2 / 3 ^/(2m). 



(9) 
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In the full crossover problem, combining the limiting val- 
ues (|6]) and we suggest the following bulk chemical 
potential in BBM valid for all negative a 
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to be used in Eq. |(7J). This is the simplest minimal 
bulk chemical potential consistent with limits © and ^j. 
Expression ([5]) also satisfies these limits and is consistent 
with the accurate GFMC calculations. We have checked 
through numerical calculation that expressions 1(5)) and 
(fTUjl agree to each other within an error of less than 1%. 
We suggest here the GL-type BBM Eqs. © and {TP]) 
including beyond mean-field corrections valid for a < 
and specially in the unitarity limit a — > — oo. 

For a > 0, from Eq. ©, in the small a' limit, the 
leading terms in energy density can be written as 
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This result, applicable to dimers, is written in terms of 
dimer variables denoted by prime, and is obtained for 
a uniform hard-sphere Bose gas (here composite bosonic 
dimers) in a perturbation calculation for small n'a' 3 [2(| . 
This leads to a bulk chemical potential 0, HH 
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in the mean-field Eq. ([7]), however, now with mass, trap, 
scattering length, particle number etc. appropriate for 
dimers. We have recovered in Eq. (| 12(1 the proper beyond 
mean-field generalization of the GP equation for small a 
[26| . In the large a limit, from Eq. (J9]), the leading term 
in the bulk chemical potential is given by [3] 



M '(n',a') =^ 2 (67r 2 ) 2/3 n' 2 / 3 /m'. 



FIG. 3: (Color online) Square of the radial and axial frequen- 
cies (a) Vp and (b) v\ in a cigar-shaped trap vs. /cfo, for BBM 
(Tlijl . GFMC, MFBCS as quoted in [p. Data for radial fre- 
quency taken from Refs. i8, 29] and those for axial frequency 
taken from Refs. [5, |30|. [3l| . 



Equations. (fT2"]) and (fTB")) can be combined to form the 
following chemical potential valid for small and large a 
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This bulk chemical potential has been constructed to sat- 
isfy Eq. fl3]) for large a' and Eq. (Jl2j) for small a'. 
(After a simple algebra it can be shown that Eq. (fl"4"|) 
satisfies limit (0) for large a' with ( « 1.) Expression 
(I14p is much simpler than that obtained directly from 
Eq. ([5]). Equations ([7]) and (fl4|) are the present GP-type 
BBM equations for the bosonic dimers including beyond 
mean-field corrections valid for a > and especially in 
the unitarity limit. For £ = 0.44, Eq. (fT4]) produces 
the following unitarity limit for composite bosonic dimers 
fjf(n', a') = nh 2 n' 2 / 3 jm! , k « 7. For fundamental bosons 
a similar relation is obtained with a different coefficient 
k . Equations ([7]) and (fT4|) with a different numerical 



coefficient £ appropriate for the study of a fundamental- 
boson superfluid has recently been suggested [27| . 

Next we subject bulk chemical potential (fl"4"]) to the 
stringent test by calculating the radial and axial frequen- 
cies Up and v z of collective oscillation in a cigar shaped 
trap, where r = (p, z) are the radial (p) and axial (z) 
coordinates. Cozzini and Stringari 0, [28| showed that 
for a power-law politropic dependence of bulk chemical 
potential on density p! cx n' T , v p and are given by 



v 



% v 2 z = 3 - l/r(r + 1), respectively, 
To test the BBM, we plot, in Fig. [3] 
(13) (a) and (b), v 2 and v 2 vs. atan(l//cFa), respectively, and 
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FIG. 4: (Color online) Density n(r) vs. r in units of oscillator 
lingth (I — ^JhjrrujS) in BCS and Bose regimes normalized to 
47T J n(r)r 2 dr = N = 2000: Unitarity \a\ — ► oo; Ideal Bose 
a — » +0; BCS a — > — 0. 

compare them with the GFMC and MFBCS calculations, 
and experimental data of Refs. [H, HH, as quoted in Refs. 
[1, Him. The end points of the present BBM plot in 
Fig. [3] are determined by the value of T — 2/3 and 1 at 
a = +oo and 0, respectively. However, as v 2 is related 
to the derivative dp! jdn! , bulk chemical potential (fT4"|) 
provides a good fit of the derivative of the // — n! curve 
to theoretical models @, [3l[ and experiment [!, [3l[ in 
addition to the p! — n' curve as can be seen from Fig. 
[3l In this figure the difference between different curves 
is small (order of 2% to 3%) and by slightly altering the 
constants in Eq. (fl"4f we can obtain a result close to 
either the GFMC or MFBCS plot. We did not do so 
and the present result is obtained only from asymptotic 
condition without any fitting. 

Finally, we solve, numerically, for a spherical harmonic 
trap, Eqs. © and (JTUJ) for fermions (a < 0) and Eqs. 
and (fT4|) for dimer bosons (a > 0) by the method of imag- 



inary time propagation after discretizing it with the semi- 
implic it Cran k-Nicholson rule. We employ length scale 
I = ^/h/muj and time scale (r = lo^ 1 ) with to the an- 
gular frequency of trap, and 2000 fermionic atoms (1000 
dimers). In numerical simulation we use I = 0.025 and 
t = 0.001. The calculated densities plotted in Fig. [4] 
show interesting behavior. The rms radius r rms in differ- 
ent regimes obey inequalities r"^~° > ri^ 00 > r"^+° 
with numerical values 4.74^m, 3.85^m, and 1.25/im, re- 
spectively. The size depends on the respective nonlin- 
earity, increasing as the nonlinearity increases. In the 
BCS limit (a — > —0) the Fermi gas extends to a greater 
distance than at unitarity due to an increased repulsion. 
The ideal dimer Bose gas with a' — has the most com- 
pact structure due to no repulsion. 

To conclude, we proposed and solved numerically 
a quasi-analytic nonlinear Eq. ([7]) for Bose to BCS 
crossover of a dilute Fermi gas with beyond mean-field 
correction so that it is valid in the unitarity region with 
divergent scattering length a. This model produces the 
known analytic behavior of the energy and bulk chem- 
ical potential of the system (dependence on scattering 
length) in the BCS (a — > —0), Bose (a — > +0), and uni- 
tarity a — > ±oo limits. For a < (Fermi regime), the 
equations are GL-type Eqs. ([7]) and (fl~0|) . and for a > 
(dimer Bose regime), they are the GP-type Eqs. ([7|) and 
(fl~4|) . The calculated radial and axial frequencies of col- 
lective breathing oscillation of a system in a cigar-shaped 
trap are found to be in good agreement with experiment 
and GFMC and MFBCS calculations. 
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